Transport in a laser irradiated thin foil 
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Three dimensional Particle-In-Cell simulations describing the interaction of a short intense laser 
pulse with thin foils are presented. It is observed that the laser generated electron current decays 
into magnetically isolated filaments. The filaments grow in scale and magnitude by magnetic re- 
connection. Two different laser wavelengths are considered. The spatial separation of the filaments 
varies for the two wavelengths. Many current filaments carry net electric currents exceeding the 
Alfven current considerably. 
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Key issues of laser-matter interaction at high intensities are the absorption of large fractions of the irradiated 
laser energy and the transport of large energy flows through a plasma. A typical application of intense laser-matter 
interaction is Fast Ignition (FI) in Inertial Confinement Fusion (ICF) jl]-[|. While FI requires large laser facilities 
high field experiments with thin foils promise to show interesting transport properties in dense plasma with moderate 
requirements for pulse energy [^,^) . Laser heated foils have been proposed for x-ray generation and harmonic emission 
HQ, for the production of well-defined dense plasma films J8|, and as targets for efficient ion acceleration |J. Hence, 
55 ■ there is both experimental and theoretical interest in a good characterization of laser- irradiated thin foils. 

We call foils thin when their thickness is larger than a few skin lengths l s = c/uj p , where uj p is the plasma frequency, 
O ■ and very much smaller than the mean free path of the electrons. This definition has some ambiguity since a reasonable 
^ ' definition depends on the physics which is investigated. The latter however is not fully known. Two dimensional (2D) 
simulations in the plane perpendicular to the laser direction have recently been reported [ p0| . In these simulations 
" 35 ' ^ ne laser has been neglected and fast particles have been injected by hand. However, it is an open question to which 
extend laser-generated spectra can be modeled. In addition, 2D geometry severely limits the available degrees of 
freedom for current transport and magnetic field evolution. 
i ^~f . Here we investigate mechanisms of self-organization in laser-generated charge flows in thin foils in full 3D where the 
laser produces the hot particle spectra. We show that merging of magnetic filaments [fi"l]l is an important process to 
. build up larger filaments. We show that this process proceeds slowly in time as soon as a balance between magnetic 
and thermal pressures in the plasma is reached. For a relatively cold plasma into which a well defined beam is 
injected it is found that the Alfven current flj;[l3| cannot be exceeded [ pj5|Jl4| . However, in this paper it is revealed 
' that laser generated charge flows can exceed a reasonably defined Alfven current limit. To address the problems 
discussed we perform two Particle-In-Cell simulations (PIC) with different parameters in 3D for a short intense laser 
pulse interacting with a sharp edged thin foil of many times over-critical plasma. Collisions will be neglected in both 
simulations. 

The simulations mentioned differ only in the laser wavelengths. The first one has A = 0.815/im and the second A = 
1.612/um. All fields depend on x, y, and z. The simulation boxes are 10/im x 10/mi x 4/im large and periodic in lateral 
• r— I ' directions. The grids have 400 x 400 x 400 cells. Electrons and ions are presented by equal quasi-particle numbers. The 
total number of quasi-particles in the simulations is 1.28 • I0 8 . The initial electron and ion temperatures are 10.0 keV 
and 1.0 keV respectively. The laser pulses propagate in z and are linearly polarized along x. After a rise time of three 
O |. optical cycles the laser intensity is kept constant. The incident laser pulses have a Gaussian envelope laterally with 
a width of 2.5/xm at full-width-half-maximum. The irradiances in both cases arc IX 2 = 6.0 • 10 18 Wcm~ 2 /im 2 . The 
foils have a marginal initial deformation to enhance absorption to generate large currents [^J. The deformations are 
' parameterized by z(x, y) = (S-exp (—(x — Xq) 2 jr 2 + (y — y ) 2 /r 2 ) where 5 = 0.4/im, xq = y$ = 5.0/Ltm, and r = 2.5/xm. 
The center of the foils are located at z = 2.2/im. The thicknesses are 0.6 /im. The background ions consist of carbon 
and are assumed to be singly charged. The plasma density in the simulations is n c = 3.33 • 10 22 cm~ 3 . The coordinate 
system used in the simulations is right handed. Hence, the positive z-axis points out of the xy-plane shown in the 
figures. In what follows we will always use the positive z-axis for reference. 

Plot (a) of Fig. |l] shows the plane y = 5.0 /zm of the cycle averaged magnetic field |B| defined in the figure caption. 
Plot (b) of the same figure gives the plane z = 2.6 /xm of |B|. Plot (c) shows a magnification of the central area of 
the foil in plot (b) at z = 2.6 /mi. The arrows plotted over the pictures indicate- the direction of the field vector of B. 
Plots (a,b,c) are for A = 0.815/mi while plots (d,e,f) show the same for A = 1.62/xm. In the center of the foil magnetic 
filaments arc observed. As is seen from plots (c,f) the magnetic field in the filaments rotates counter-clockwise when 
viewed along the positive z-axis. The peak magnetic field strength which is obtained in the filaments is about 8700 T 
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for A = 0.815^m and 5500 T for A = 1.63/im. The white arrows in plots (b,e) show that the central magnetic filaments 
are surrounded by a magnetic fields of larger scale. The topology of this field is such that it rotates clockwise in front 
of the foil where z < 2.2/^m holds and counter-clockwise at the rear of it for z > 2.2 /im. This is a consequence of 
electrons escaping into the vacuum both at the front and rear of the foil in opposite directions. The electrons escape 
from the center and return along the front and rear foil surfaces. The total electric current obtained by integrating 
over the whole lateral foil extension disappears on both the front and rear sides of the foil (see Fig. ||). 

Plots (a,b,c) of Fig. show the evolution of the magnetic filaments for A = 1.63/um. Plot (a) indicates the early 
stage of the filament evolution. Filaments are very small initially. There have been efforts to describe this early stage 
with the help of the Weibel theory (ljjl^l ■ Later the small filaments reconnect H] to form larger scales as is seen from 
plot (b). Plot (c) shows the saturation scale. Reconnection in the present context can be understood qualitatively 
with the help of Eq. (|^) derived later in the paper which describes the total force on a finite plasma volume V. In 
case V contains a single filament for which j x B has a large degree of asymmetry a net force is obtained. When 
two filaments approach each other the magnetic fields of the latter tend to cancel in the overlapping region leading 
to further asymmetry in j x B. Hence, the two filaments reconnect to form a larger one. The magnetic field topology 
is changing in the course of the process since initially closed field lines in different filaments end up as a single new 
filament with again closed field lines. The energy required is provided by the laser driver. 

Plot (a) of Fig. H shows the cycle averaged current filaments obtained in the foil for A = 0.815/xm. The simulation 
reveals roughly twenty of them where only the dark dots visible in the center of plot (a) have been counted for a 
filament. Plot (b) of the same figure gives a magnified view of the central current filaments for better illustration. 
The peak current density in the filaments is \j z \ = 1.6 • 10 17 A/m 2 . Plots (c,d) show the same for A = 1.62/im. The 
peak current density in the filaments is now \ j z \ = 7.6 • 10 16 A/m 2 . It is evident that the filament separation and scale 
are larger for A = 1.62/im while the current density is approximately half the short wavelength value. Current density 
and filament scale are correlated as is implied by Eq. (|l]) derived later in this paper. 

Plots (a,b,c) of Fig. |^ compare current density, total current, and Alfven current in the central region of the foil 
for A = 1.62/xm. The white rectangles indicate the area over which the current density in (a) has been integrated 
to obtain the total current given in (b). The size of the rectangle is such that it matches approximately the cross 
sectional area of a single filament. Plot (b) has been obtained moving the rectangle over plot (a) point-wise. The 
Alfven current shown in (c) has been calculated with the help of Ia = 17.0 PjkA, where (5 = v/c and 7 = 1/ \Jl — /3 2 . 
The quantity (3 has been obtained from the averaged velocity obtained from the fastest particles that constitute the 
net current in a particular space-time point. This means that c/3(x,i) = ff ast d 3 pv z /(x, p, t)/ J d 3 p /(x, p, t). Next, 
/3(x, t) | z =2.6/im is averaged over the cross sectional area indicated by the white rectangle shown in plot (c). It is seen 
that the Alfven current limit thus obtained is exceeded substantially by the net current present in the highlighted 
current filament. Wc note that there are large ambiguities about how to calculate the Alfven current limit in a hot 
laser plasma since no unique particle beam can be identified in the particle spectrum. We note further that in our 
simulations fast electrons do not contribute significantly to the net charge flow. The latter is predominantly due to 
slow electrons. As a consequence the averaged f3 is small. 

Plot (a) of Fig. H shows the electric and magnetic field energy distributions for both simulations obtained from 
the cycle averaged fields E and B by integrating the field energy densities laterally over the whole foil extension. 
In the vacuum the electric field energy dominates. In the foil the magnetic field energy dominates while the cycle 
averaged electric field disappears. To show that the total current and return current balance while locally they differ 
substantially plot (b) of the same figure shows the current in both foils obtained by integrating j z laterally as in plot 
(a). The peak values of the total electric currents are approximately / = 5.0 ■ 10 5 A for A = 0.815/^m and I = 2.5 • 10 5 A 
for A = 1.62/^m. For comparison, the total currents limited by the Alfven current are 7.6 ■ 10 4 A for A = 0.815/im and 
4.1 • 10 4 A for A = 1.63/xm. 

To derive a relation between the magnetic force acting on a filament and the thermal pressure inside the foil we 
neglect the electric field E since it disappears there. In addition, we are only interested in slowly varying time scales. 
Under these assumptions we find making use of the Vlasov equation 



Summation over repeated indices is implied. The quantity P l} denotes the pressure tensor which is obtained directly 
from the simulations. The integral form of Eq. ([j]) is given by 





(1) 




(2) 
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According to Eq. (g) the force on the plasma volume V is obtained by the pressure acting on the surface surrounding 
the volume. From the simulations we find that |6|jP z j| <C |i9jP X j|, |<9jP y j| holds. This is consistent with the observation 
that |j x B| x . y 3> |j x B| z holds inside the foil. Figure || shows magnitude and direction of 9jPj obtained from the 
simulations. Multiplying the magnetic field obtained from plot (f) of Fig. [l] with the current density j z in plot (d) of 
Fig. m and comparing with Fig. ^ it is seen that the magnetic force on the current is approximately balanced by the 
kinetic pressure of the electrons for the quasi-steady state shown in the present paper. 

Fig. ^ shows the electron spectrum obtained for A = 1.62/im. We find that current and return current contain 
fast electron populations. The fast electrons in the return current are fast electrons that have been reflected at the 
back surface due to the electric field present there. Comparing plots (a,b) shows that the fast electrons in the return 
current try to escape from the filaments. Peak electron energies obtained for A = 1.63/rai are of the order of 1.0 MeV 
as is seen from plots (a,b) while the electron temperature of the hot electrons found within the inner rectangle (see 
plot (c) of Fig. R and plot (a) of Fig. ^) as well as within the inner and outer rectangles (see plot (d) of Fig. ^ and 
plot (a) of Fig. [|) is roughly 300.0 kcV. To understand these low values we note that the foil represents a sharp edged 
plasma with n/n c = 80.0 for A = 1.63/im. Hence, the laser radiation decays rapidly over distances of l s ~ 0.03/im. 
As a consequence the effective fields in the plasma capable of accelerating electrons become very small. 

Transport properties in thick foils cannot be directly predicted on the basis of those found in thin foils. The main 
reason is that the electron distribution function in thin foils tends to become symmetric in momentum space which is 
not so for extended plasma slabs. As is found in this paper, aa approximate force balance between thermal pressure 
and magnetic force can be established in thin foils which is capable of slowing down further filament reconnection. It 
needs to be investigated two which degree the latter is true in thick plasma slabs. 

In conclusion, we have shown with the help of 3D PIC simulations that large charge flows generated by intense 
laser radiation in thin foils decay into current filaments. The latter originate from the critical surface and extend into 
the bulk plasma. Their separation and magnitude is correlated with the laser wavelength. The current filaments are 
surrounded by magnetic filaments. The evolution of the magnetic filaments is mainly governed by filament merging. 
In our simulations this process slows down when an approximate balance between thermal pressure and the magnetic 
force acting on the filaments is reached. The net charge flow in a single filament can exceed a reasonable definition 
of the Alfven current substantially. For the short wavelength case studied in this paper the net current in a single 
filament exceeds the Alfven current up to 8-fold. However, a reasonable Alfven current definition turns out to be 
ambiguous in the context of hot laser plasmas. 

This work has been supported by the European Commission through the TMR network SILASI, contract No. 
ERBFMRX-CT96-0043. Use of the Supercomputing facilities at ZIB (Konrad Zuse, Berlin, Germany) and NIC (John 
von Neumann Institute, Jiilich, Germany) has been made. 
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FIG. 1. Cycle averaged B. Plot (a) shows the plane y — 5.0 pm, plots (b,c) the planes z = 2.6 /im of the cycle averaged 
magnetic field B = (Bl + By + B\ ) 6 for A = 0.81pm and t = 54 fs. Plots (d,e,f) show the same for A = 1. 62pm and t = 130 fs. 
The arrows indicate the direction of the cycle averaged magnetic field B. The parameter is Bo = 8.76 ■ 10 2 Vs/m 2 . 

FIG. 2. filament merging in B. The plots show the planes z = 2.6/tm of the cycle averaged magnetic field 
B = {Bl + Bl+Blf 5 for A = 1.62pm and t = 30fs (a), t = 40fs (b), and t = 108 fs (c). The parameter is B = 8.76-10 2 Vs/m 2 . 

FIG. 3. Current density j z . Plots (a,b) show the planes z — 2.6 pm of the cycle averaged electron current density j z for 
A = 0.81/im and t = 54 fs. Plots (c,d) show the same for A = 1.62/iin and t = 130 fs. The arrows indicate the direction of the 
cycle averaged magnetic field B. The parameter is jo = 1.66 • 10 16 A/m 2 . 

FIG. 4. Alfven current 7a- The plots show the planes z = 2.6/im of the cycle averaged current density j z (a), the total 
current I obtained from integrating j z over the white square (b), and the Alfven current 7a- The actual current exceeds the 
Alfven current 4-fold. The parameters are jo = 1.66 ■ 10 16 A/m 2 , A = 1.62/iin, and t = 130fs. 

FIG. 5. Field energy and total current 7. Plot (a) shows the electromagnetic field energies per micron in the foils calculated 
from the time averaged fields. Dotted lines correspond to A = 0.81/im and t = 54fs. Solid ones belong to A = 1.62/im and 
t — 130fs. Bold lines give the magnetic field energy for both cases which has been multiplied by a factor of five. Plot (b) shows 
the electron current 7 Z . Dotted lines correspond to A = 0.81/im and t — 54fs. Solid ones belong to A = 1.62/im and t = 130fs. 

FIG. 6. Pressure gradient. Cycle averaged thermal force density |9j J-\j|. The colorbar is in units j B = 1.46 • 10 19 N/m 3 for 
comparison with Figs. |l| and ^. The white arrows show the force direction. The parameters are A = 1.62/im and t = 130fs. 

FIG. 7. Electron spectrum. Plot (a) shows the zp z -plane of the phase space located inside the inner rectangle indicated in 
plot (a) of Fig. |^. Plot (b) shows the same for particles located between the inner and outer rectangles in plot (a) of Fig. ^. 
Plot (c) shows the energy spectrum obtained form the particles located in the inner rectangle while plot (d) shows the one 
obtained from electrons between the inner and outer rectangle. The parameters are A = 1.62pm and t = 130fs. 
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